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Abstract 

This paper is devoted to investigate gravitational collapse of dust 
in metric f{R) gravity. We take FRW metric for the interior region 
while the Schwarzchild spacetime is considered for the exterior region 
of a star. The junction conditions have been derived to match interior 
and exterior spacetimes. The assumption of constant scalar curvature 
is used to find a solution of field equations. Gravitational mass is 
found by using the junction conditions. It is concluded that the con¬ 
stant curvature term f{Ro) plays the role of the cosmological constant 
involved in the field equations of general relativity. 
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1 Introduction 


The most interesting topic in the gravitational physics today is the expansion 
of our universe. The support to this argument comes from different sources 
such as Supernovae la (SNIa) experiments [1], cosmic microwave background 
fluctuations |2] and X-ray experiments |3]. All these observations indicate 
that our universe is expanding with an accelerated rate. The phenomenon of 
dark matter and dark energy is another topic of discussion |1] . Einstein gave 
the concept of dark energy in 1917 by introducing a small positive cosmo¬ 
logical constant in the held equations. But later on, he rejected it. However, 
it is now believed that the cosmological constant may be a suitable candi¬ 
date for dark energy. Higher dimensional theories [S] such as M-theory or 
string theory may also be helpful to explain this cosmic expansion. Another 
justihcation comes from modihcation of general theory (GR) involving some 
inverse curvature terms [6]. However, modihed gravity with inverse curvature 
terms seems to be unstable and may not pass solar system tests [7] . This dis¬ 
crepancy can be addressed by including higher derivative terms. Moreover, 
the viability can be achieved by considering squared curvature terms [8]. It 
has been suggested that the current expansion may be justihed if we add 
some suitable powers of curvature in the usual Einstein-Hilbert action |9]. 
Thus it seems interesting to study the universe in the context of alternative 
or modihed theories of gravity. 

Among various modihcation, f{R) gravity is a possible candidate which 
gives a natural gravitational alternative to dark energy [10]. This theory may 
provide an easy unihcation of early time inhation and late time acceleration. 
The cosmic acceleration can be explained by introducing the term 1/i? at 
small curvatures. It was Buchdahl m who introduced f{R) gravity using 
non-linear Lagrangians. The f{R) theory of gravity seems most suitable due 
to its cosmologically important f{R) models. These models consist of higher 
order curvature terms as functions of Ricci scalar R. Some viable f{R) grav¬ 
ity models |9] have been suggested which show the unihcation of late-time 
acceleration and early-time inhation. It is now expected that dark matter 
problem can be addressed using viable f{R) models. In recent years, many 
authors have shown keen interest to investigate this theory in diherent con¬ 
text m-m- Some detailed reviews are available to better understand the 
theory [18]. Multamaki and Vilja [19] investigated static spherically symmet¬ 
ric vacuum solutions in f{R) theory. They established that the held equations 
in f{R) gravity gave the Schwarzschild de Sitter solution. Exact spherically 
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symmetric interior solutions in metric f{R) gravity have been studied by 
Shojai and Shojai [20] • Hollenstein and Lobo |2I] analyzed static spherically 
symmetric solutions in f{R) gravity coupled to non-linear electrodynamics. 
f{R) gravity at one-loop level in de Sitter universe has been investigated 
by Cognola et al. [22]. Cylindrical symmetry has also been widely used to 
investigate f{R) gravity in different contexts [23] • Recently developed /(T) 
gravity is another alternative theory which is the generalization of teleparal¬ 
lel gravity. This theory also seems interesting as it may explain the cosmic 
acceleration without involving the dark energy. A considerable amount of 
work has been done in this theory so far [23]- 

Gravitational collapse is an interesting and important issue in GR. It is 
involved in the structure formation of the universe causing the existence of 
galaxies, stars and planets. The singularity theorem suggests that the occur¬ 
rence of spacetime singularity is a general feature of any cosmological model 
under some reasonable conditions. So the solutions with singularities can 
be produced by the gravitational collapse of non-singular and asymptotically 
flat initial data [23]. The classification of spacetime singularities is based on 
two facts whether they can be observed or not. If a spacetime singularity 
is locally observable then it is termed as naked. A black hole is a space- 
time singularity which can not be observed. Penrose [23] proposed a cosmic 
censorship conjecture that the singularities appearing in the gravitational col¬ 
lapse are always covered by an event horizon. Formation of compact stellar 
objects like neutron stars and white dwarf is the result of gravitational col¬ 
lapse. Spherical symmetry has been extensively used to study gravitational 
collapse. Dust gravitational collapse was first explored by Oppenheimer and 
Snyder [27]. Markovic and Shapiro [2E] extended their work by consider¬ 
ing positive cosmological constant. Many researchers [22] have investigated 
gravitational collapse by considering interior and exterior regions. 

In the recent years, many authors have shown keen interest to explore 
gravitational collapse in alternative theories of gravity [30]- [33]. It has been 
shown that farther one goes from GR, there is a greater chance of having a 
naked singularity [3D] • Ghosh and Maharaj [33] obtained a condition for the 
occurrence of a naked singularity in the collapse of null dust in higher dimen¬ 
sional f{R) gravity. Openheimer-Snyder collapse in Brans-Dicke theory has 
been discussed by Scheel [33]. In a recent paper [37], Rudra and Debnath 
discussed gravitational collapse in Vaidya spacetime for Galileon theory of 
gravity. Sharif and Abbas [33] studied the dynamics of shearfree dissipative 
gravitational collapse in f{G) gravity. Spherically symmetric perfect fluid 
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gravitational collapse has been discussed in metric f{R) gravity by Sharif 
and Kausar [35] • Cembranos et ah jlU] analyzed a general f{R) model with 
uniformly collapsing cloud of self-gravitating dust particles. 

In this paper, we are focussed to discuss the gravitational collapse with 
dust case in f{R) gravity. We take Friedmann-Robertson-Walker (FRW) 
spacetime in the interior region and Schwarzschild metric in the exterior 
region. The paper is organized as follows: Section 2 is used to give general 
formalism about junction conditions between interior and exterior regions. 
We introduce held equations in f{R) gravity and solve them using FRW 
metric for dust case in section 3. In Section 4, we hnd the apparent horizons 
and discuss the role of constant curvature term. Finally, we summarize the 
results in the last section. 


2 General Formalism 

In this section, we give junction conditions at the surface of a collapsing dust 
sphere. For this purpose, a 4D spherically symmetric spacetime is divided by 
a time-like 3D hypersurface S into two regions namely interior and exterior 
regions. Interior and exterior regions are denoted by V~ and V~^ respectively. 
Interior region represented by FRW spacetime is given by 

= dt‘^ — a^{t)dr‘^ — {t)b‘^ {r)[d6‘^ + Sin‘^6d(j)^], (1) 

where a(t) is cosmic scale factor and 

{ sinr, when k=l, 

r, when k=0, 

sinhr, when k=-l. 

For exterior region V~^, we consider the Schwarzschild spacetime 

dsl = (1 - —)dT^ -+ Sin^ed(j)^], ( 2 ) 

^ 1- -fT 

where M is an arbitrary constant. Using Israel junction conditions, we con¬ 
sider that hrst and second fundamental forms for interior and exterior space- 
times are same. These conditions are given as: 

1. The continuity of hrst fundamental form over S provides 

{dsD^ = {ds^_h = {ds%, (3) 
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2. The continuity of second fundamental form over S yields 


[Kij] =^tj - Kij, {i, j = 0, 2, 3), 


where the extrinsic curvature tensor Kij is dehned as 




± 




± 


r 


dx^ dx'i 


de'^de^ de'^ de^ 


{a,ii,u = 0,1,2,3). 


(4) 

(5) 


Here e* and correspond to the coordinates on S and respectively. 
The christoffel symbols are calculated using the interior and exterior 
spacetimes and are the components of outward unit normal to S in the 
coordinates x^. Using interior and exterior spacetimes, the equations of 
hypersurface S are written as 


h_(r, t) = r — Vs = 0, 
h+{R,T) = R-R^{T) = 0, 


( 6 ) 

(7) 


where rs is an arbitrary constant. Using these equations, interior and exterior 
metrics given in Eq.(IT]) and(l2]) take the form 


{dslh = 


1 - 


(ds^)s = dt^ — a'^{t)b‘^{rj:)[d6‘^ + Sin‘^6d(j)‘^], 
2M 1 


R, 


^ — Ml 

i?E 



dT^ - Rl 

\dT J \ 



d9^ + Sin^ed^^ 


( 8 ) 


(9) 


Here we assume 1 — ^ > 0 so that T remains time-like 

coordinate. Using junction condition (|3]), we get 

2 


1 - 


2M 


R, 


dRs 


1 _ flV dT 

ht . 


dT = dt. 


( 10 ) 

( 11 ) 


Now using Eqs. ([6]) and ([7]), the outward unit normals to the interior and 
exterior spacetimes are given by 


= (0,a(t),0,0), 


( 12 ) 


5 

















( 13 ) 


n+ = (-i?s,T,0,0), 

while the components of extrinsic curvature turn out to be 
Kqo = 0, K 22 = csc^OK^^ = (a66')s, 

3Mi?2t M{R-2M)f^^ 


K^o = 


RT-TR 


R{R - 2M) 

K+ = csc^eK+ = {f{R - 2M))s. 

Here dot and prime denote differentiation with respect to t and r respectively. 
By applying continuity condition on extrinsic curvatures, we obtain 


( 14 ) 

( 15 ) 

( 16 ) 


^00 ~ 0 ; 


— K~ 
^22 — ^22- 


(17) 


Now using Eqs. (11411171) along with (HOI) and (Dll. the junction conditions 
take the form 


{b% = 0, M 


+ a(Rb^ — abb''^'^ 


(18) 


It is mentioned here that equations (ITOl) . (ITTI) and (|T8|l forms necessary and 
sufficient conditions to match the interior and exterior regions smoothly. 


3 f{R) Gravity and Field Equations 

The metric tensor has a key role in GR. One of the main features of GR 
is the dependence of Levi-Givita connection on the metric tensor. However, 
the connection does not remain the Levi-Givita connection if we allow the 
torsion in the theory. Gonsequently the dependence of connection on metric 
tensor vanishes. This is the main idea behind different approaches of f{R) 
theories of gravity. 

We get metric version of f{R) gravity if the connection is the Levi-Givita 
connection. In this approach, the variation of action is done with respect to 
the metric tensor only. The action for f{R) gravity is 

SfiR) = j + Lm)d\ (19) 

where f{R) is a general function of the Ricci scalar and Lm is the mat¬ 
ter Lagrangian. It would be worthwhile to mention here that the standard 
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Einstein-Hilbert action can be achieved when f{R) = R. Varying this action 
with respect to the metric tensor yields the modihed held equations 

F{R)R^, - ^f{R)g,u - + gf,.DF{R) = kTJZ- (20) 

Here F{R) = df{R)/dR, n is the coupling constant, T™ is the standard 
energy-momentum tensor and 




( 21 ) 


with is the covariant derivative. We can write the held equations in an 
alternative form which is familiar with GR held equations 

- \a„uR = + f™, (22) 


where = TJ^/F{R) and the energy-momentum tensor for gravitational 
huid is 


F{R) 


-g,u\fiR)-RF{R) ) +FiR) 


;q;/3 


9a}i9pu 9^v9a^ 


( 23 ) 


It can be seen from Eq. (1231) that energy-momentum tensor for gravitational 
huid contributes matter part from geometric origin. This approach is 
interesting as it may provide all the matter components required to investi¬ 
gate the dark universe. Thus it is hoped that f{R) theory of gravity may 
give fruitful results to understand the phenomenon of expansion of universe. 
When we contract Eq. fEU]) . it follows that 


F{R)R - 2f{R) + 3DF{R) = kT™, (24) 


where T™ is the trace of energy-momentum tensor. Here we are interested in 
pressureless gravitational collapse. For dust, the energy-momentum tensor 
is given as 

= pUf^u^, (25) 

where p is the matter density and the four velocity vector satishes the 
equation Using this equation along with held equations fl20|) . we 

get three independent diherential equations for the interior spacetime 


3a 

a 


1 

F 




(26) 
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( 27 ) 


a 

a 


a , d, 9 h" 1 
_ + 2{-)^ - 2— = — 
a a?b F 


(1 1 


f aF ^ 

- + 2-h F 

2 a 


1 

F 


/ + 2^ + f 

2 a 


(28) 


The solution of these highly non-linear differential equations does not seem 
to be possible straightforwardly. However, we can try to find a solution 
using the assumption of constant scalar curvature, i.e. R = Rq. Using this 
assumption, left side of Eq. fl2T|) becomes constant which leads to a constant 
energy density, say p = po. Thus the held equations fl261l28|) take the form 


3d 

a 


F{Ro 


Kpo 


- + 2A^-2^ = 

a a a^o 


a 

a 


h' 

a?h a?F 


,a 2 F 


f{Ro 

f{Ro) 

2F{R^y 

1 f{Ro) 


2F{Ro 


Manipulating these equations, we obtain 

1 


,a ,a 2 1 


2“ + (“) + 2h2 

a a FF 


2F{Ro 

Integrating hrst equation from (|T8ll . it follows that 

b’ = X, 


Mo + /(Bo)]' 


(29) 

(30) 

(31) 

(32) 

(33) 


where X in an arbitrary integration function of r. Thus Eq. fl3^ takes the 
form 


2- + {-Y + 

a a 


FF 2F(i?o) 

Integrating this equation with respect to t, we get 

+ fiRo) 


. 2 — 1 ^ 2m 

“ F ^ ^ 


6F{R 


'O) L 


(34) 


(35) 


where m = m{r) is an arbitrary function and is related to the mass of the 
collapsing system 


m{r) 


upoFF 

QF{Ro)' 
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(36) 

























Using gravitational units, i.e. k 
takes the form 

m{r) 


Stt, the mass of the collapsing system 


4:7rpoa%^ 

3F(i?o) ■ 


It is mentioned here that mass of the system must be positive because neg¬ 
ative mass is not acceptable physically. Using Eq. (l35ll and second junction 
condition in Eq. ffTSD . we get 


M = m — 


a^&^[87rpo -F f{Ro)] 

12F[Rq) 


(38) 


Now we calculate the total energy M{r,t) at a time t for the interior hyper¬ 
surface of radius r using the mass function |5T] 


Using Eq. 


, the mass function turns out to be 

a^6^[87rpo + f{Ro)] 


M{r, t) = m{r) 


12F{Ro 


(39) 


(40) 


where m(r) denotes the energy due to constant matter density in Eq. fl37)) . 
Now we hnd the solution with X(r) = 1 using Eq. fl33|l . In this case the closed 
form solution turns out to be 


ab 


— 12mF(i?o) 

Sttpo + f{Ro) 


2 

sinhs a{r, t), 


(41) 


where 


a{r, t) = 


3[87rpo + f{Ro)] 


8F{Rq 


[u 


t]. 


(42) 


Here we assume 87rpo + f{Ro) < 0 to have a realistic solution and ts{r) is 
an arbitary function of r. It is clear that t = tg is the time formation of 
singularity for a perticular shell at some distance r. In the limiting case 
when /(i?o) ^ the above solution takes the form 


lim ab 

f{Ro)^-8iTpo 



( 43 ) 


which correspond to the well known Tolman-Bondi solution |42] . 
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4 Apparent Horizons 

We obtain the apparent horizon when the bonndary of two trapped spheres 
is formed. In this section, we hnd snch bonndary of two trapped spheres 
whose ontward normals are nnll. For the interior spacetime ([1]), this is given 
as 


{ah) = a^b — 6'^ = 0. 

Using Eq. fl35D in this eqnation, we get 


(44) 


—[Svrpo + f{Ro))]a%^ + 6a6 - 12m = 0. 
-c irtoj 


(45) 


The solntions of this eqnation for ab yield the apparent horizons. For f{Ro)) = 
—87ipo, it becomes the Schwarzschild horizon, i.e., ab = 2m. When m = 0, 
it yields a de-Sitter horizon, i.e. 



(46) 


The case 3m < 



) Is^-ds to two horizons. 



(47) 


and 



(48) 


where the snbscripts c and bh represent cosmological and black hole horizons 
respectively and 



(49) 


If we take m = 0, the equations fjTTl) and fj48ll reduce to 



( 50 ) 
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It is mentioned here that the results can be generalized when m 7 ^ 0 and 
Svrpo + fiRo) 7^ 0 


For the case when 3m = sj both horizons coincide, i.e., 

{ab)c = {ab)hh = 


-2F{Ro 


(51) 


Sttpo + fiRo)' 

Thus the range of the cosmological horizon and the black hole horizon be¬ 


comes 


0 < iab)bh < 


-2F(i?o 


< iab)c < 


- 6 F(i?o 


(52) 


Svrpo + fiRo) \j 8npo -|- fiRo) 

There does not exist any apparent horizon in the case 3m > sir 

The formation time of the apparent horizon can be calculated using Eqs. 
flTTl) and fjT5l) and is given by 


t'n. to 


— 8 F(i?o 


3(87rpo + fiRo)) 


sinh 


-1 


iab)„ 

2 m 


- 1 


u = 1,2. 


(53) 


In the limiting case when fiRo) —STipo, the result corresponds to Tolman- 
Bondi solution 


From Eq. fl5^ . it follows that 


where 


a„(r,f) = 


4 

tah 3 ^* 

(54) 

ictb)n ^2 

^ = cosh On, 

2m 

(55) 

-3(8^p„+ /(«,)) 

SFiRo) ^ ^ ^ ' 

(56) 


It is clear from Eq. fl53p that both the black hole horizon and the cosmological 
horizon form earlier than the singularity t = tg- Now since Eq. fl5^ yields 
ti < t 2 and Eq. (l5^ implies that iab)bh < iab)c- This is an indication that 
cosmological horizon forms earlier than the black hole horizon. 
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5 Concluding Remarks 


This paper is devoted to discuss gravitational collapse in f{R) gravity. For 
this purpose, we consider the metric approach of this theory to study the field 
equations. It is observed that the held equations fl26l) - fl28l) are complicated 
and highly non-linear. Thus it seems difficult to solve them analytically 
without any assumption. The assumption of constant curvature has been 
used to hnd a solution. 

We investigate the gravitational collapse of dust by considering FRW 
spacetime as the interior region while for the exterior region we take Schwarzschild 
metric. The junction conditions have been derived between interior and exte¬ 
rior spacetimes. We get two physical apparent horizons namely cosmological 
horizon and black hole horizon. It is found that formation time for black hole 
horizon is more as compared to cosmological horizon. Moreover, both hori¬ 
zons are formed earlier than singularity. This indicates that the singularity is 
covered, i.e., black hole, which shows that f{R) gravity supports cosmic cen¬ 
sorship conjecture. From dynamical equation 0351) . the rate of gravitational 
collapse is 


The acceleration should be negative for the collapsing process which is pos¬ 
sible when 


ab < 


6 mF(i?o) 
Sttpo + /(^o) 


1 

3 


(58) 


Thus Eq. (I57|) indicates that /(i?o) slows down the collapsing process when 
/(i?o) + Sttpo < 0 . Further, the presence of f{Ro) causes two physical hori¬ 
zons. One is the black hole horizon and the other is cosmological horizon. 
It also influences the time difference between the formation of the apparent 
horizon and singularity. It is concluded that f{Ro) affects the process of 
collapse and hence it limits the size of the black hole. It would be worth¬ 
while to mention hare that the term f{Ro) plays the same role as that of the 
cosmological constant in GR held equations and our results agree with |44) . 


12 







Acknowledgement 


MFS is thankful to National University of Computer and Emerging Sciences 
(NUCES) Lahore Campus, for funding the PhD programme. The authors 
are also grateful to the anonymous reviewer for valuable comments and sug¬ 
gestions to improve the paper 


References 

[1] Perlmutter, S. et al.: Astrophys. J. 517(1999)565; Riess, A.G. et 
al.; Astron. J. 116(1998)1009; Astier, P. et al.: Astron. Astrophys. 
447(2006)31. 

[2] Seprgel, D. N. et al.: Astrophys. J. Suppl. Ser. 148(2003)175; Bennett, 
C. L. et ah: Astrophys. J. 148(2003)1. 

[3] Allen, S. W. et al.: Mon. Not. R. Astron. Soc. 353(2004)457. 

[4] Nojiri, S. and Odintsov, S.D.: Int. J. Geom. Meth. Mod. Phys. 
4(2007)115; Turner, M.S., Huterer, D.: J. Phys. Soc. Jap. 
76(2007)111015; Frieman, J., Turner, M. and Huterer, D.: Ann. Rev. 
Astron. Astrophys. 46(2008)385; Li, M., Li, X.D., Wang, S. and Wang, 
Y.: Commun. Theor. Phys. 56(2011)525; Copeland, E.J., Sami, M. 
and Tsujikawa, S.: Int. J. Mod. Phys. 015(2006)1753; Sahni, V. and 
Starobinsky, A.: Int. J. Mod. Phys. D9 (2000) 373; Sahni, V.: Lect. 
Notes. Phys. 653 (2004) 141; Carroll, S.M.: Living Rev. Rel. 4(2001)1; 
Weinberg, D.H.: New. Astron. Rev. 49(2005)337; Straumann, N.: Mod. 
Phys. Lett. A21 (2006) 1083. 

[5] Nojiri, S. and Odintsov, S.D.: Phys. Lett. B576(2003)5; Guenther, 
U., Zhuk, A., Bezerra, V. and Romero, C.: Class. Quant. Grav. 
22(2005)3135. 

[6] Carroll, S.M., Duvvuri, V., Trodden, M. and Turner, M.: Phys. Rev. 
D70(2004)043528; Capozziello, S., Carloni, S. and Troisi, A.: Int. J. 
Mod. Phys. D12(2003)1969. 

[7] Chiba, T., Phys. Lett. B575(2003)l; Soussa, M.E. and Woodard, R.P.: 
Gen. Rel. Grav. 36(2004)855. 


13 



[8] Nojiri, S. and Odintsov, S.D.: Phys. Rev. D68(2003) 123512; Abdalla, 
E., Nojiri, S. and Odintsov, S.D.: Class. Quant. Grav. 22(2005)L35. 

[9] Nojiri, S. and Odintsov, S.D.: Phys. Rev. Nojiri, S. and Odintsov, S.D.; 
Problems of Modern Theoretical Physics, A Volume in honour of Prof. 
Buchbinder, I.L. in the occasion of his 60th birthday, p.266-285, (TSPU 
Publishing, Tomsk), arXiv:0807.0685, 

[10] Nojiri, S. and Odintsov, S.D.: Int. J. Geom. Meth. Mod. Phys. 
115(2007)4. 

[11] Buchdahl, H.A.: Mon. Not. Roy. Astr. Soc. 150(1970)1. 

[12] Sharif, M. and Shamir, M.F.: Glass. Quantum Grav. 26(2009)235020; 
Sharif, M. and Shamir, M.F.: Mod. Phys. Lett. A25(2010)1281; Sharif, 
M. and Shamir, M.F.: Gen. Relativ. Gravit. 42(2010)2643; Sharif, M. 
and Zubair, M.: Adv. High Energy Phys. 2013(2013)790967; Sharif, M. 
and Kausar, H.R.; JGAP 07(2011)022; Sharif, M. and Kausar, H.R.; 
Int. J. Mod. Phys. 020(2011)2239. 

[13] Bamba, K., Nojiri, S., Odintsov, S.D. and Saez-Gomez, D.: Phys. Lett. 
B730(2014)136; Bamba, K., Makarenko, A.N., Myagky, A.N., Nojiri, 
S. and Odintsov, S.D.: JGAP 01(2014)008; Bamba, K., Nojiri, S. and 
Odintsov, S.D.: Phys. Lett. B698(2011)451; Nojiri, S. and Odintsov, 
S.D.: Phys. Rept. 505(2011)59; Elizalde, E., Nojiri, S., Odintsov, S.D. 
and Saez-Gomez, D.: Fur. Phys. J. 070(2010)351; Bamba, K., Geng, 
G., Nojiri, S. and Odintsov, S.D.: Mod. Phys. Lett. A25(2010)900; 
Nojiri, S. and Odintsov, S.D.: Phys. Rev. D68(2003)123512. 

[14] Gapozziello, S. and Vignolo, S.: Int. J. Geom. Meth. Mod. Phys. 
8(2011)167; Gapozziello, S., Darabi, F. and Vernieri, D.: Mod. Phys. 
Lett. A26(2011)65; Gapozziello, S., Laurentis, M.D., Odintsov, S.D. 
and Stabile, A.: Phys. Rev. D83(2011)064004; Gapozziello, S., Lauren¬ 
tis, M.D., Nojiri, S. and Odintsov, S.D.: Gen. Rel. Grav. 41(2009)2313; 
Gapozziello, S. and Vignolo, S.: Glass. Quantum Grav. 26(2009)175013. 

[15] Paul, B.G., Debnath, P.S. and Ghose, S.: Phys. Rev. D79(2009)083534. 

[16] Shamir, M.F., Jhangeer, A. and Bhatti, A.A.: Ghin. Phys. Lett. 
29(8)(2012)080402. 


14 


[17] Jamil, M., Mahomed, F.M and Momeni, D.; Phys. Lett. 
B702(2011)315; Momeni, D., Raza, M and Myrzakulov, R.: Eur. Phys. 
J. Plus 129(2014)30; Farooq, M. U., Jamil, M., Momeni, D. and Myrza¬ 
kulov, R.: Can. J. Phys. 91(2013)703; Hendi, S.H. and Momeni, D.; 
Eur. Phys. J. 071(2011)1823. 

[18] Sotiriou, T.P. and Faraoni, V.: Rev. Mod. Phys. 82(2010)451; Felice, 
A.D and Tsujikawa, S.: Living Rev. Rel. 13(2010)3; Nojiri, S. and 
Odintsov, S.D.; Phys. Rept. 505(2011)59; Bamba, K., Capozziello, S., 
Nojiri, S. and Odintsov, S.D.: Astrophys. Space Sci. 342(2012)155. 

[19] Multamaki, T. and Vilja, L: Phys. Rev. D74(2006)064022. 

[20] Shojai, A. and Shojai, F.: Gen. Relativ. Gravit. 44(2011)211. 

[21] Hollenstein, L. and Lobo, F.S.N.; Phys. Rev. D78(2008)124007. 

[22] Cognola, G., Elizalde, E., Noiiri, S., Odintsov, S.D. and Zerbini, S.: 
JCAP 0502(2005)010. 

[23] Rodrigues, M.E., Houndjo, M.J.S., Momeni, D. and Myrzakulov, R.: 
Can. J. Phy. 92(2)(2014)173; Azadi, A., Momeni, D. and Nouri-Zonoz, 
M.: Phys. Lett. B670(2008)210; Momeni, D. and Gholizade, H.; Int. J. 
Mod. Phys. D18(2009)l; Sharif, M. and Arif, S.: Astrophys. Space Sci. 
342(2012)237. 

[24] Jamil, M., Momeni, D. and Myrzakulov, R.: Eur. Phys. J. 
072(2012)1959; ibid. 072(2012)2075; ibid. 072(2012)2137; ibid. 
073(2013)2267; Jamil, M., Yesmakhanova, K., Momeni, D. and Myrza¬ 
kulov, R.: Cent. Eur. J. Phys. 10(2012)1065; Houndjo, M.J.S., Mo¬ 
meni, D. and Myrzakulov, R.: Int. J. Mod. Phys. D21 (2012) 1250093; 
Yang, R.J.: Europhys. Lett. 93(2011)60001; Wei, H., Ma, X.P. and 
Qi, H.Y.: Phys. Lett. B703(2011)74; Wu, P.X. and Yu, H.W.: 
Eur. Phys. J. 071(2011)1552; Wu, P.X. and Yu, H.W.: Phys. Lett. 
B703(2011)223; Bamba, K., Geng, C.Q., Lee, C.C and Luo, L.W.: 
JCAP 1101(2011)021. 

[25] Penrose, R.; Phys. Rev. Lett. 14(1965)57; Hawking, S.W.: Proc. R. 
Soc. London A300(1967)187. 

[26] Penrose, R.: Riv. Nuovo Cimento 1(1969)252. 


15 



[27] Oppenheimer, J.R. and Snyder, H.: Phys. Rev. 56(1939)455. 

[28] Markovic, D. and Shapiro, S.L.: Phys. Rev. D61(2000)084029. 

[29] Misner, C.W. and Sharp, D.: Phys. Rev. D136(1964)b571; Ghosh, 
S.G. and Deshkar, D.W.: Int. J. Mod. Phys. 012(2003)317; Sharif, 
M. and Ahmad, Z.: J. Korean Phys. Society 52(2008)980; Sharif, M. 
and Ahmad, Z.: Int. J. Mod. Phys. A23(2008)181; Sharif, M. and 
Abbas, G.: J. Phys. Soc. Jpn. 80(2011)104002; Sharif, M. and Ab¬ 
bas, G.: Gen. Relativ. Gravit. 44(2012)2353; Sharif, M. and Abbas, 
G.: Astrophys. Space Sci. 327(2010)285; Sharif, M. and Abbas, G.: 
J. Phys. Soc. Jpn. 82(2013a)034006; Sharif, M. and Abbas, G.: Ghi- 
nese Phys. B22(2013)030401; Sharif, M. and Abbas, G.: Eur. Phys. J. 
Plus 28(2013)10; Abbas G. and Ramzan, R.M.: Ghinese Phys. Lett. 
30(2013)100403; Debnath, U., Nath, S. and Ghakraborty, S.: Mon. Not. 

R. Astron. Soc. 369(2006)1961; Debnath, U., Nath, S. and Ghakraborty, 

S. : Gen. Relativ. Grav. 37(2005)215. 

[30] Debnath, U., Ghakraborty, N. G., Ghakraborty, S.: Gen. Rel. Grav. 
40(2008)749;Debnath, U., Ghakraborty, S., Barrow, J. D. : Gen. Rel. 
Grav. 36(2004)231; Banerjee, A.,Debnath, U., Ghakraborty, S.dnt. J. 
Mod. Phys. D12(2003)1255; Rudra, P., Biswas, R. and Debnath, U.: 
Astrophys Space Sci. 335(2011)505; Debnath, U., Rudra, P. and Biswas, 

R. : Astrophys Space Sci. 339(2012)135; Rudra, P., Biswas, R. and Deb¬ 
nath, U.: Astrophys Space Sci. 342(2011)557. 

[31] Maeda, H.: Phys. Rev. D73(2006) 104004. 

[32] Jhingan, S. and Ghosh, S.G.: Phys. Rev. D81(2010)024010; Ghosh, 

S. G. and Jhingan, S.: Phys. Rev. D82(2010)024017. 

[33] Joshi, P.S. and Singh, T.P. : Phys.Rev. D51 (1995)6778; Singh, T.P. and 
Joshi, P.S.: Glass. Quant. Grav. 13(1996)559; Joshi, P.S. and Dwivedi, 
I.H.; Gommun. Math. Phys. 146(1992)333; Joshi, P.S. and Dwivedi, 
I.H.; Lett. Math. Phys. 27(1993)235. 

[34] Path, K.D. and Thool, U.S.; Int. J. Mod. Phys. D14(2005)873. 

[35] Ghosh, S.G. and Maharaj, S.D.: Phys. Rev. D85(2012)124064. 

[36] Scheel, M.A. et ah: Phys. Rev. D51(1995)4236. 


16 



[37] Rudra, P. and Debnath, U. :arXiv:1402.0350vl. 

[38] Sharif, M. and Abbas, G.: J. Phys. Soc. Jpn. 82(2013)034006. 

[39] Sharif, M. and Kansar, H.R.: Astrophys. Space Sci. 331(2011)281. 

[40] Cembranosa, J.A.R, Cruz-Dombrizb, A.D.L. and Nnnez B.M.: JCAP 
04(2012)021. 

[41] Misner, C.W. and Sharp, D.; Phys. Rev. 136(1964)b571. 

[42] Eardley, D.M. and Smarr, L.: Phys. Rev. D19(1979)2239. 

[43] Hayward, S.A., Shiromizn, T. and Nakao, K.: Phys. Rev. 

D49(1994)5080. 

[44] Sharif, M. and Ahmad, Z.: Mod. Phys. Lett. A22(2007)1493; Sharif, 
M. and Ahmad, Z.: World App. Sci. J. 16(2010)1516. 


17 


